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Genetic variants affecting the neural processing 
of human facial expressions: evidence using a 
genome-wide functional imaging approach 

AA Brown^'2'^\ J Jensen^ YS Nikolova^ S Djurovic^ '^ I Agartz^ ^ A Server^ RE Ferrell^ SB Manuck^ M MattingsdaP'^°, 
I Melle^'^, AR Hariri^, A Frigessi^ and OA Andreassen^'^ 

Human faces present crucial visual information for social interaction. Specialized brain regions are involved in the perception of 
faces, with the fusiform face area (FFA) a key neuronal substrate. Face processing is genetically controlled, but by which specific 
genes is unknown. A genome-wide approach identified common single nucleotide polymorphisms (SNPs) associated with areas 
of increased brain activity in response to affective facial expressions, measured with functional magnetic resonance imaging. 
SNPs in 20 genetic regions were linked with neural responses to negative facial expressions in a Norwegian sample (n = 246), 
which included patients with mental illness. Three genetic regions were linked with FFA activation in a further discovery 
experiment using positive facial expressions and involving many of the same individuals (/? = 284). Two of these three regions 
showed significant association with right FFA activation to negative facial expressions in an independent North American 
replication sample of healthy Caucasians (n = 85, 3q26.31, P= 0.004; 20p12.3, P= 0.045). The activation patterns were 
particularly striking for the SNP in 3q26.31, which lies in a gene TMEM212] only the FFA was activated. The specialized function 
of this brain region suggests that TMEM212 could contribute to the innate architecture of face processing. 
Translational Psychiatry (20^2) 2, e143; doi:1 0.1 038/tp.201 2.67; published online 24 July 2012 



Introduction 

Face processing is a crucial cognitive ability. Rapid recogni- 
tion of faces identifies individuals and informs us about basic 
properties such as age, gender and race, as well as complex 
processes such as direction of attention, mood and intent. 
Thus, the ability to extract such information from even a 
momentary viewing of a face is critically important for normal 
social interactions. 

Several lines of evidence indicate that face processing is 
unique and qualitatively different from processing of other 
visual stimuli. ^'^ Studies of neurological disorders suggest a 
dedicated circuitry exists in the human brain for face process- 
ing. Converging evidence now suggests that a network of 
multiple brain regions, in particular the right fusiform gyrus, 
supports face perception. ^"^ 

Similar specialized processing of faces has been reported 
in macaques,^'^ and it has been suggested that specialized 
pathways for face processing were critical for the survival of 
our primate ancestors.^ This implies that face processing is 
under some genetic control, and findings in infants show that 
face recognition is an innate human capability.''^''''' In fact, 
using behavioral tasks several twin studies have detected a 



high degree of heritability in face processing''^ ''^ (one study 
found that correlation of one particular score of face recog- 
nition for monozygotic twins (0.70) was more than double that 
of dizygotic twins correlation (0.29)^^). However, the specific 
genes involved in face processing are unknown, as is the 
effects of such genes on the neural networks supporting this 
specialized behavior. 

The aim of the current study was to identify common genetic 
variants associated with neural activation during face proces- 
sing. Using genome-wide genotype microarray technology 
together with functional magnetic resonance imaging (fMRI) 
in an ethnically homogenous sample of healthy controls and 
psychiatric patients from NonA/ay,'"'^ and an independent replica- 
tion sample of Caucasians from North America,''^ we discovered 
two genomic regions associated with neural activation during 
face processing. 

Materials and methods 

The Norwegian TOP sample consisted of 246 individuals after 
quality control on the genotype and imaging data. The fMRI 
protocol consisted of a widely utilized emotional faces para- 
digm, ''^"''^ involving the perceptual processing (matching) of 
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emotional facial expressions (experimental condition) or geo- 
metric shapes (control condition). Subjects completed two 
sessions, one with all negative (angry or afraid) and another 
all positive (happy) expressions, in a counterbalanced order. 
We will refer to these experiments as the negative and positive 
faces experiments. No one was excluded based on perfor- 
mance as it was a very easy task and subjects practiced it 
to demonstrate understanding before entering the scanner. 
In the negative faces experiment, the age was 34 ±10 
(mean ± s.d.) years and there were 1 30 men and 1 1 6 women. 
There were 138 participants with mental disorder, whereas 
108 were healthy controls. For the positive faces paradigm, 
we included 284 individuals, age 34 ± 9 years; 1 57 male, 1 27 
female, 160 with mental disorder and 124 healthy controls. 
Single nucleotide polymorphism (SNP) genotyping was 
performed using the Affymetrix Gene Chip Genome Wide 
SNP 6.0 array (Affymetrix Inc., Santa Clara, CA, USA, 
http://www.affymetrix.com). Individuals with discrepancies 
between reported and genotype-inferred gender and those 
whose ancestry differed from the majority of the sample were 
removed. Mitochondrial and X-chromosome SNPs were also 
removed, and we excluded individuals and SNPs with >3% 
missing data, and SNPs with a minor allele frequency of 
<0.05. 

The fMRI data from both tasks were processed using the 
SPM2 software''^ to produce contrast maps. For both negative 
and positive faces experiments, we focused on the 20,000 
voxels with most evidence of differential activation to reduce 
multiple testing. A further prescreening stage rejected voxels 
in clusters of <50. In the negative faces experiment, every 
voxel was tested for association with every SNP using an 
additive model of genetic effect. For each voxel, the P-value 
for the most significant SNP was adjusted for multiple testing 
across the SNPs using the Sidak correction. Multiple testing 
across the voxels was accounted for by controlling the false 
discovery rate (FDR) under the assumption of positive 
dependencies using a common adjustment method. ^° We 
validated our associations in the positive faces experiment by 
testing each voxel for association with 35 SNPs suggested by 
the negative faces experiment; multiple testing across voxels 
was controlled for each SNP separately.^° In both analyses, 
voxels were only included in the results if they lay within 
clusters of > 10. 

The North American replication sample consisted of 85 
healthy middle-aged Caucasian subjects, 41 males, 44 
females; age 45 ± 7 years). All fMRI data preprocessing and 
single-subject analyses were completed using SPM2.''^ 
These participants were genotyped using the lllumina 610- 
Quad BeadChip (lllumina Inc., San Diego, CA, USA). Perfect 
proxies were available on this array in two of the genetic 
regions under study, while the third region had no proxy and 
so was separately genotyped using fluorescence polarization 
methods. Single-subject contrast images were entered into 
separate regression analyses with the rsl 2485367, rs6038686 
and rs6081495 genotypes, including gender and age as 
nuisance covariates; these analyzes were conducted in 
SPM8^^ (http://www.fil.ion.ucl.ac.uk/spm). 

More details, including a full description of the subjects, and 
the processing and analysis of the data can be found in 
Supplementary material available online. 



Results 

Faces task in Norwegian sample. We carried out genome- 
wide analyses using the Affymetrix 6.0 Array on the subjects 
who had successfully completed the fMRI protocol. Clusters 
of voxels exhibiting significant differential activation in the 
direct comparisons of the experimental and control condi- 
tions were used. The 20,000 most significant voxels were 
identified and any clusters of <50 voxels were discarded, 
in order to choose groups of voxels able to contain diffuse 
associations. After this procedure, 19,742 voxels remained 
for analysis from the negative faces scan, and 19,864 from 
the positive faces scan. For both the negative and positive 
faces, voxels exhibiting a significant main effects of task (that 
is, experimental > control condition) were principally located 
in ventral visual cortices as well as frontal and subcortical 
regions. There is highly significant evidence of differential 
activation for these voxels (Bonferroni-corrected family-wise 
error (FWE) rate for the voxels in the negative faces is 
1.2x1 0~^). The exact location of the voxels selected for both 
scans is provided in Supplementary Figure SI online. 

Genetic variants suggested by negative faces 
experiment. Initially, fMRI data from the scan with nega- 
tive facial expressions were analyzed as these stimuli have 
been more commonly used than positive expressions in the 
larger neuroimaging literature.^'' Each of the voxels as 
chosen above was tested for an association with each SNP 
(a total of 549,640 SNPs, meaning approximately 10''° tests) 
using an additive model of genetic effect. The P-values were 
adjusted to account for multiple testing: first using the Sidak 
correction''^ to produce P-values, adjusted for multiple testing 
across the SNPs. Then the Benjamini Hochberg procedure^° 
was applied to these adjusted P-values to control the FDR 
(the procedure controls the expected proportion of false 
positive results), accounting for multiple testing across both 
SNPs and voxels. With the exception of the preselection of 
task-related voxels described above, this is a similar strategy 
to one previously presented applied in a genome-wide asso- 
ciation study of structural imaging phenotypes.^^ In that 
paper, identical procedures for controlling for multiple testing 
across SNPs and phenotypes (vertices as opposed to 
voxels) were followed, with the exception that we applied a 
more conservative procedure to control for multiple testing 
across SNPs (ignoring linkage disequilibrium structure and 
correcting for every SNP rather than calculating an 'effective 
number of tests'). Similarly, a set of phenotype-genotype 
associations is presented, with the appropriate estimate of 
the FDR. More details can be found in Supplementary material, 
including Supplementary Figure S2 available online, a visual 
representation of the correction procedure. 

Evidence that the activation of 455 voxels had a genetic 
component was found, with an associated FDR of 0.5. This is 
a highly liberal threshold, but as this first part of the study was 
intended to contribute to the initial discovery potential SNPs 
and significance will be assessed in the replication experi- 
ments, a higher proportion of false positives was tolerated to 
avoid discarding interesting results. 

Only clusters of at least 10 voxels were selected. This is a 
smaller threshold than that used for preselection of the 20,000 



Translational Psychiatry 



Genetic variants and human facial expressions 

AA Brown et al 



voxels (which insisted clusters had at least 50 voxels), as here 
we simply wished to discard spurious findings. There were 1 0 
such clusters, totaling 226 voxels, located in the cuneus and 
occipital gyrus, both known to be involved in visual proces- 
sing, and bilateral fusiform gyrus, which is implicated as 
central to face processing specifically.^ These clusters were 
associated with 37 SNPs located in 20 genomic regions 
<40kb. These 20 genetic regions and the associated brain 
regions are listed in Supplementary Table S1 available online. 
Supplementary Figures S3-S5 available online show Man- 
hattan and Q-Q plots for three voxel phenotypes, those most 
significantly associated with SNPs rs1 2485367, rs1 6992973 
and rs2208796, these are SNPs we will develop further. There 
is little evidence that any of these SNPs are not in Hardy 
Weinberg Equilibrium (for rs1 2485367, P= 1.00; rs1 6992973, 
P=0.14; and rs2208796, P=0.13). We see little evidence 
of genomic inflation for these phenotypes, for example, 
A. = 0.997 for the phenotype associated with rs1 2485367. 
Supplementary Figures S6-S8 available online shows in 
more detail the P-values for these phenotypes in the genetic 
regions where they show most significant associations. In 
Supplementary Figure S6 available online, we see that the 
association relates to a single signal, as no other SNPs in high 
linkage disequilibrium with rs1 2485367 were available in this 
dataset. This increases the chance that the association could 
be a false positive induced by genotyping errors; to rule this 
out, we have explored the population genetics of this SNP in 
comparison with the North American sample and the CEU 
sample from the 1 000 genomes pilot project^^ for evidence of 
biased ascertainment. This information is contained in the 
Supplementary Material. We also note that in Supplementary 
Figure S6 available online, we see that a number of SNPs 
closest to rs1 2485367 show P-values between 10~^ and 
10""^ despite being in only weak linkage disequilibrium with 
rs1 2485367, some evidence that they may be weakly tagging 
the causal variant. 

We plot genotype against peak phenotype for these three 
SNPs in Supplementary Figure S9 available online. All three 
SNPs seem to show a consistent trend across genotype; there 
is little evidence these associations are induced by a few 
extreme outliers. Finally, Supplementary Figure S1 0 available 
online shows Manhattan and Q-Q plots for the peak voxel 
associated with each of the 20 genetic regions listed in 



Supplementary Table S1 online. Phenotypes from certain 
brain regions, in particular the cuneus, show more evidence of 
genomic inflation. This could be because these regions show 
greater variation in brain activation, possibly breaking the 
normality assumptions of the model. 

Validation of results in positive faces experiment. The 

37 most significant SNPs from the negative faces scan were 
tested in the positive faces scan. This is not an independent 
replication in a new population, as the experiments were 
performed on mostly the same people during the same 
scanning occasion performing a similar task and thus the 
P-values reported are biased towards significance. However, 
recruitment into the study continued after data for the nega- 
tive faces scan were collated and analyzed and so genotype 
data on 57 new individuals were available for the positive 
faces scan. Two SNPs were discarded as their minor allele 
frequency was <0.05 in this new population, leaving 35 SNPs. 
Again using an additive model of genetic effect, associations 
between SNPs in seven genetic regions and neural activation 
were identified in 26 clusters of >10 voxels, located in the 
cuneus, lingual and superior occipital gyri as well as in the 
bilateral fusiform gyrus (the associated voxel-wise FDR 
was 0.05, the associations with regions of the fusiform 
gyrus, a structure particularly involved in face processing, are 
listed in Table 1, all associations discovered are listed in 
Supplementary Table S2 online). 

The neural regions associated with the seven genetic 
regions across both the negative and positive faces data are 
all known to be involved in visual processing. Some of these 
regions, such as the primary occipital gyrus, are involved with 
general visual processing. Others overlap with the so-called 
occipital face area which, together with the fusiform gyrus, 
is considered necessary for normal face processing. For 
the fusiform gyrus, the right hemisphere is known to be 
preferentially involved in face processing.^ The results in the 
negative faces scan are consistent with this, a greater number 
of associations were found with fusiform gyrus voxels in the 
right hemisphere than the left hemisphere (36 voxels in the 
right fusiform gyrus were associated with genetic variants, as 
opposed to 19 in the left fusiform gyrus). The pattern of 
activation in the Norwegian sample is particularly striking in 
the case of rsl 2485367; significant associations with the 



Table 1 SNPs associated with the fusiform gyrus for the positive faces tasl< 



SNP 


Gene annotation 


Location 
(minor aileie frequency) 


Hemisphere 


Pea/f coordinate 


Cluster size 


Unadjusted 
peal< significance 


rs2208796 


SLC24A3 


20:19048148 (0.19) 


Left 


-22 -62 -2 


52 


3.79 X 10"^ 


rsl 2485367 


TMEM212 


3:173114366 (0.069) 


Right 


38 -50-10 


47 


1.00 X 10"^ 






Left 


-36 -52 -12 


31 


7.66 X 10"^ 








Left 


-44 -44-12 


12 


2.94 X 10"^ 


rsl 6992973 


RP4-764O22-001 


20:6957158 (0.058) 


Left 


-34 -62 -10 


21 


6.87 X 10"^ 






Left 


-30 -44 -24 


14 


7.40 X 10"^ 








Right 


32 -68 -10 


14 


5.52 X 10"^ 



Abbreviation: SNP, single nucleotide polymorphism. 

Three genetic regions from the negative faces task were also associated with activation with the fusiform gyrus in the positive faces task, located in seven clusters of 
voxels. We present the names and locations of the most significant SNP, the hemisphere within which the association is seen and the most significant nominal P-value 
for association. These associations were further investigated in an independent sample. 
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fusiform gyrus, only in the right hemisphere in the negative 
faces scan (1 6 voxels), and more pronounced in the right than 
left hemisphere for the positive faces scan were discovered 
(linked to 47 voxels in the right hemisphere and 43 in the left). 
This SNP was not associated with any other brain regions in 
either experiment. In addition, the peak activation in the right 
hemisphere to negative faces is in almost exactly the same 
location as the activation associated with this SNP in the 
positive faces scan (see Figure 1). In the positive faces 
experiment, the rs1 6992973 variant was also associated 
with bilateral fusiform area activation (though more pro- 
nounced in the left hemisphere); the remaining genetic 
region was associated only with activation in the left 
fusiform gyrus. 

Our initial analysis ignored disease status, as analyses of 
pilot data found no trend effects on activation. However, as 
others have found associations between regional brain 
activation and mental illness, it was important to investigate 
whether the present results could be caused by confounding 
with this factor. We found no significant associations between 
disease status and the activation of the 20,000 preselected 
voxels in either experiment. We have also investigated how 
much controlling for disease status, using indicator variables 
coding for schizophrenia, bipolar and other psychosis, would 
affect the results. Supplementary Figure S1 1 available online 
shows the P-values declared significant in both experiments 
against the same P-values calculated when diagnosis is 
included in the model (for the negative faces experiment, all 
associations listed in Supplementary Table S1 available 
online were reanalyzed, this time regressing on SNP dose 
and disease status; for the positive faces experiment, all 
associations listed in Supplementary Table S2 available 
online were analyzed similarly). We see little change; for 
both experiments, results were actually slightly improved 
when diagnosis was included in the model, with a small 
majority of associations more significant (in the negative faces 
experiment, 207 associations were more significant includ- 
ing diagnosis, 197 were more significant ignoring it; for the 
positive faces experiment, these numbers were 627 and 615, 
respectively). 




Figure 1 Pattern of activation for rs1 2485367 in TMEM212. We present coronal 
sections -52), showing tliis single nucleotide polymorphism to be associated 
with activation in fusiform gyrus; unilaterally in the negative faces paradigm 
displayed on the left, and bilaterally in the positive faces paradigm (panel to right). 
Activated voxels are circled and colored red. The nominal threshold for significance 
is 2. 1 X 1 0"^ for the negative faces paradigm and 3.2 x 1 0"^ for the positive faces 
paradigm; after adjusting for multiple testing corresponds to a false discovery rate of 
0.5 and 0.05. 



We have also repeated the full analysis performed on the 
positive faces experiment, but again, controlling for diagnosis. 
Similarly, we found that the results were largely unchanged, 
indeed slightly improved. Considering diagnosis makes very 
little difference to the conclusions we draw and would not 
affect the choice of SNPs to be investigated further in the 
replication experiment described below. This agrees with 
previous work, which found preserved function of the fusiform 
face area (FFA) in schizophrenia.^^ These results are covered 
in more detail in the Supplementary material online. 

Replicating fusiform gyrus activations in Nortli American 
sample. Associations across the negative and positive faces 
datasets were identified. However, in the first negative faces 
dataset, we expect half of the associations to be false posi- 
tives, and in the second positive faces dataset, the estimates 
of significance are difficult to interpret as the SNPs were 
selected using related data. Therefore, and in line with recent 
guidelines for replicating genotype-phenotype associations,^^ 
a subset of the significant SNPs was investigated in an 
independent sample of Caucasians from North America. The 
replication experiments were designed to concentrate on 
those SNPs most likely to be involved in face processing. 
The different levels of complexity of the stimuli could cause 
activations; viewing faces could also inspire higher levels of 
attention than geometric figures as could the emotional 
content. However, there is a substantial body of literature 
pointing to a specialized role for the FFA in face processings^ 
and thus we concentrate our focus on SNPs linked to this 
brain region, in particular in the positive faces paradigm 
(one SNP of particular interest, rs1 2485367, was associated 
with FFA activation in both negative and positive faces 
paradigms). We use the literature to argue that activations in 
this region are related to face processing rather than any 
other cognitive subtractions. For SNPs associated with other 
regions, a replicated finding would not have this body of work 
to aid in the interpretation. 

The SNPs that were identified in the negative faces dataset 
and associated with FFA activation in the positive faces 
dataset were selected (there were five such SNPs, located in 
three regions). The analyses on activations in the fusiform 
gyrus used an anatomical region of interest approach, using 
standard fMRI analyses, which are impractical in the genome- 
wide association study context due to the large number of 
tests. The NCI-NHGRI Working Group guidelines state that 
replication experiments and analyses should be as similar as 
possible, though this is chiefly in order to remove study 
differences as a reason for failure to replicate. 

The replication sample consisted of healthy middle-aged 
volunteers who completed a similar fMRI task''^''^ from which 
the NoHA/egian paradigms were derived. These participants 
were genotyped using the lllumina platform. The SNPs 
upstream of SLC24A3 (20p1 1 .23), had a proxy on the lllumina 
platform (rs6081495). Another SNP, rs1 6992973 near RP4- 
764022-001 (20p12.3), had a proxy in rs6038686. All proxies 
were in complete linkage disequilibrium (PF and D equal 1 in 
the European ancestry population of the HapMap project^^). 
The final SNP, rs1 2485367 in TMEM212 (3q26.31), did not 
have a proxy on the lllumina platform and was individually 
genotyped. 
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A significant association between rs1 2485367 (TMEM212) 
and right fusiform gyrus activation was obtained (x=20, 
y=-68, z=-^4, t=4.77, P=0.004, FWE-corrected, /ce=92; 
see also Figure 2). As in the original experiments, participants 
carrying the minor (G) allele showed higher fusiform gyrus 
activation relative to C homozygotes. Significant association 
between rs6038686 (RP4-764O22-001) and right fusiform 
gyrus activation (x= 26, y= -52, z= -1 4; T= 4.07, P= 0.043, 
FWE-corrected, /ce=80; see Figure 3) was also found; such 
that participants carrying the minor (G) allele showed rela- 
tively heightened activation. Again, this matches the result 
from the Norwegian sample. As in the rs1 2485367 analyses, 
the right fusiform gyrus activation remained significant after 
FWE correction over the search volume of an anatomical 
fusiform gyrus region of interest. Finally, rs6081495 did not 
show a significant association with fusiform gyrus activation. 

Discussion 

The present findings represent a demonstration of an 
association between common gene variation and neural 
function associated with face processing. Validation of some 
of the initial results using a parallel experiment and replication 
of the most promising candidates in an independent popula- 
tion was used. It is also important to note that the original 
associations with fusiform gyrus activation in the NonA/egian 
discovery sample were identified when searching across all 



neural regions sensitive to the task and not through the narrow 
window of an anatomical region of interest approach. 

There is little knowledge about the function of the genes that 
were identified, but there are some indications that TMEM212 
is expressed in the brain. Affymetrix provides a test of whether 
a transcript is present or absent within a sample,^^ and 
TMEM212 is present in five out of seven samples of the whole 
brain tissue collated in the Gene Enrichment Profiler^° (http:// 
xavierlab2.mgh.harvard.edu/EnrichmentProfiler/index.html), 
and in one out of two samples from the occipital lobe. 

Interestingly, both genes replicated in the North American 
sample lie within regions for which associations have been 
proposed to disorders involving impairment of social function- 
ji^g 31-35 g^j^^ association to diagnosis was found, 

and the replication sample excludes any individuals with a 
current DSM-IV (Diagnostic and Statistical Manual for Mental 
Disorders) Axis I disorder, we believe these genes to be 
involved in regulating fundamental brain functions related to 
face processing unrelated to psychiatric disorders. Further, as 
these variants were identified and validated using faces 
expressing negative and positive emotions, these pathways 
may be broadly important for face processing independent of 
emotional expression. 

The current findings are in line with recent evidence of high 
heritability of the ability to process face information,''^'''^ which 
suggests a strong biological control of this important human 
function. Previous studies indicate that the ability for face 




Figure 2 Replication of an association between TMEM212 and fusiform gyrus activation in a North American Caucasian sample. Significant association between 
rs1 2485367 and right fusiform gyrus activation is presented in coronal, sagittal and axial overlays on a high-resolution single-subject T1 -weighted structural image. Participants 
carrying the minor (G) allele exhibited relatively heightened activation in comparison to those homozygous for the C allele. Activation parameters: cluster location in MNI 
(Montreal Neurological Institute) space: x=20, y= -68, z= -14, 7=4.77, P= 0.004, FWE (family-wise error)-corrected, /c£=92. 
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Figure 3 Replication of an association between RP4-764O22-001 and fusiform face area response in a North American Caucasian sample. Significant association 
between rs6038686 and right posterior fusiform response is presented in coronal, sagittal and axial overlays on a high-resolution single-subject T1 -weighted structural image. 
Participants carrying the minor (G) allele exhibited relatively heightened activation in comparison to those homozygous for the A allele. Activation parameters: cluster location 
in MNI (Montreal Neurological Institute) space: x=26, y= -52, z= -14; 7=4.07, P= 0.043, /ce=80. 



recognition has little overlap with other cognitive abilities, ''^ ''^ 
which suggests that the face processing ability is specific. This 
is supported by the present findings of specific activation 
patterns related to face processing associated with distinct 
gene variants, indicating a specific genetic basis. Particularly 
TMEM212 seems to be specifically involved in the innate 
architecture of face processing, as it was associated with 
activation in the FFA only. 

This supports the position that distinct cognitive and 
neuronal mental processes may have specific genetic deter- 
minants, as has been seen with language disability.^^ Thus, 
the genes identified in the current study would be interesting 
candidate genes to investigate for association with develop- 
mental prosopagnosia,^^ a disorder characterized by impair- 
ments in face recognition that are unaccompanied by brain 
lesions; multiple cases of this disorder have been seen in 
the same family. ^^'^^ Investigations into pedigrees has 
suggested that forms of congenital prosopagnosia can have 
a dominant, autosomal mode of inheritance. Although we 
have looked for common variation, these genes could be 
investigated to find rare variants causing pathological 
results in such families. However, the development of 
such distinct mental processes probably depends on gene- 
environment interactions. 

To conclude, specific genes that may be important for the 
modulation of basic neurobiological processes, which under- 
lie human social interactions, have been identified. In addition, 
we have demonstrated how the synthesis of fMRI and 



genome-wide data can be used to discover and replicate 
links between common genomic variation and behaviorally 
relevant brain function. 
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